The mechanism shaping the logistic growth of mutation proportion in epidemics at population scale.

Virus evolution is a common process of pathogen adaption to host population and environment. Frequently, a small but important fraction of virus mutations are reported to contribute to higher risks of host infection, which is one of the major determinants of infectious diseases outbreaks at population scale. The key mutations contributing to transmission advantage of a genetic variant often grow and reach fixation rapidly. Based on classic epidemiology theories of disease transmission, we proposed a mechanistic explanation of the process that between-host transmission advantage may shape the observed logistic curve of the mutation proportion in population. The logistic growth of mutation is further generalized by incorporating time-varying selective pressure to account for impacts of external factors on pathogen adaptiveness. The proposed model is implemented in real-world data of COVID-19 to capture the emerging trends and changing dynamics of the B.1.1.7 strains of SARS-CoV-2 in England. The model characterizes and establishes the underlying theoretical mechanism that shapes the logistic growth of mutation in population.

SARS-CoV-2 lineage dynamics in England from September to November 2021: high diversity of Delta sub-lineages and increased transmissibility of AY.4.2.

BACKGROUND: Since the emergence of SARS-CoV-2, evolutionary pressure has driven large increases in the transmissibility of the virus. However, with increasing levels of immunity through vaccination and natural infection the evolutionary pressure will switch towards immune escape. Genomic surveillance in regions of high immunity is crucial in detecting emerging variants that can more successfully navigate the immune landscape. METHODS: We present phylogenetic relationships and lineage dynamics within England (a country with high levels of immunity), as inferred from a random community sample of individuals who provided a self-administered throat and nose swab for rt-PCR testing as part of the REal-time Assessment of Community Transmission-1 (REACT-1) study. During round 14 (9 September-27 September 2021) and 15 (19 October-5 November 2021) lineages were determined for 1322 positive individuals, with 27.1% of those which reported their symptom status reporting no symptoms in the previous month. RESULTS: We identified 44 unique lineages, all of which were Delta or Delta sub-lineages, and found a reduction in their mutation rate over the study period. The proportion of the Delta sub-lineage AY.4.2 was increasing, with a reproduction number 15% (95% CI 8-23%) greater than the most prevalent lineage, AY.4. Further, AY.4.2 was less associated with the most predictive COVID-19 symptoms (p = 0.029) and had a reduced mutation rate (p = 0.050). Both AY.4.2 and AY.4 were found to be geographically clustered in September but this was no longer the case by late October/early November, with only the lineage AY.6 exhibiting clustering towards the South of England. CONCLUSIONS: As SARS-CoV-2 moves towards endemicity and new variants emerge, genomic data obtained from random community samples can augment routine surveillance data without the potential biases introduced due to higher sampling rates of symptomatic individuals.

The non-pharmaceutical interventions may affect the advantage in transmission of mutated variants during epidemics: A conceptual model for COVID-19.

As the COVID-19 pandemic continues, genetic mutations in SARS-CoV-2 emerge, and some of them are found more contagious than the previously identified strains, acting as the major mechanism for many large-scale epidemics. The transmission advantage of mutated variants is widely believed as an innate biological feature that is difficult to be altered by artificial factors. In this study, we explore how non-pharmaceutical interventions (NPI) may affect transmission advantage. A two-strain compartmental epidemic model is proposed and simulated to investigate the biological mechanism of the relationships among different NPIs, the changes in transmissibility of each strain and transmission advantage. Although the NPIs are effective in flattening the epidemic curve, we demonstrate that NPIs probably lead to a decline in transmission advantage, which is likely to occur if the NPIs become intensive. Our findings uncover the mechanistic relationship between NPIs and transmission advantage dynamically, and highlight the important role of NPIs not only in controlling the intensity of epidemics but also in slowing or even containing the growth of the proportion of variants.

Estimating the transmission advantage of the D614G mutant strain of SARS-CoV-2, December 2019 to June 2020.

IntroductionThe SARS-CoV-2 lineages carrying the amino acid change D614G have become the dominant variants in the global COVID-19 pandemic. By June 2021, all the emerging variants of concern carried the D614G mutation. The rapid spread of the G614 mutant suggests that it may have a transmission advantage over the D614 wildtype.AimOur objective was to estimate the transmission advantage of D614G by integrating phylogenetic and epidemiological analysis.MethodsWe assume that the mutation D614G was the only site of interest which characterised the two cocirculating virus strains by June 2020, but their differential transmissibility might be attributable to a combination of D614G and other mutations. We define the fitness of G614 as the ratio of the basic reproduction number of the strain with G614 to the strain with D614 and applied an epidemiological framework for fitness inference to analyse SARS-CoV-2 surveillance and sequence data.ResultsUsing this framework, we estimated that the G614 mutant is 31% (95% credible interval: 28-34) more transmissible than the D614 wildtype. Therefore, interventions that were previously effective in containing or mitigating the D614 wildtype (e.g. in China, Vietnam and Thailand) may be less effective against the G614 mutant.ConclusionOur framework can be readily integrated into current SARS-CoV-2 surveillance to monitor the emergence and fitness of mutant strains such that pandemic surveillance, disease control and development of treatment and vaccines can be adjusted dynamically.

Real-time quantification of the transmission advantage associated with a single mutation in pathogen genomes: a case study on the D614G substitution of SARS-CoV-2.

BACKGROUND: The COVID-19 pandemic poses serious threats to global health, and the emerging mutation in SARS-CoV-2 genomes, e.g., the D614G substitution, is one of the major challenges of disease control. Characterizing the role of the mutation activities is of importance to understand how the evolution of pathogen shapes the epidemiological outcomes at population scale. METHODS: We developed a statistical framework to reconstruct variant-specific reproduction numbers and estimate transmission advantage associated with the mutation activities marked by single substitution empirically. Using likelihood-based approach, the model is exemplified with the COVID-19 surveillance data from January 1 to June 30, 2020 in California, USA. We explore the potential of this framework to generate early warning signals for detecting transmission advantage on a real-time basis. RESULTS: The modelling framework in this study links together the mutation activity at molecular scale and COVID-19 transmissibility at population scale. We find a significant transmission advantage of COVID-19 associated with the D614G substitution, which increases the infectivity by 54% (95%CI: 36, 72). For the early alarming potentials, the analytical framework is demonstrated to detect this transmission advantage, before the mutation reaches dominance, on a real-time basis. CONCLUSIONS: We reported an evidence of transmission advantage associated with D614G substitution, and highlighted the real-time estimating potentials of modelling framework.

The emergence and ongoing convergent evolution of the SARS-CoV-2 N501Y lineages.

The independent emergence late in 2020 of the B.1.1.7, B.1.351, and P.1 lineages of SARS-CoV-2 prompted renewed concerns about the evolutionary capacity of this virus to overcome public health interventions and rising population immunity. Here, by examining patterns of synonymous and non-synonymous mutations that have accumulated in SARS-CoV-2 genomes since the pandemic began, we find that the emergence of these three "501Y lineages" coincided with a major global shift in the selective forces acting on various SARS-CoV-2 genes. Following their emergence, the adaptive evolution of 501Y lineage viruses has involved repeated selectively favored convergent mutations at 35 genome sites, mutations we refer to as the 501Y meta-signature. The ongoing convergence of viruses in many other lineages on this meta-signature suggests that it includes multiple mutation combinations capable of promoting the persistence of diverse SARS-CoV-2 lineages in the face of mounting host immune recognition.

Quantification of the spread of SARS-CoV-2 variant B.1.1.7 in Switzerland.

BACKGROUND: In December 2020, the United Kingdom (UK) reported a SARS-CoV-2 Variant of Concern (VoC) which is now named B.1.1.7. Based on initial data from the UK and later data from other countries, this variant was estimated to have a transmission fitness advantage of around 40-80 % (Volz et al., 2021; Leung et al., 2021; Davies et al., 2021). AIM: This study aims to estimate the transmission fitness advantage and the effective reproductive number of B.1.1.7 through time based on data from Switzerland. METHODS: We generated whole genome sequences from 11.8 % of all confirmed SARS-CoV-2 cases in Switzerland between 14 December 2020 and 11 March 2021. Based on these data, we determine the daily frequency of the B.1.1.7 variant and quantify the variant's transmission fitness advantage on a national and a regional scale. RESULTS: We estimate B.1.1.7 had a transmission fitness advantage of 43-52 % compared to the other variants circulating in Switzerland during the study period. Further, we estimate B.1.1.7 had a reproductive number above 1 from 01 January 2021 until the end of the study period, compared to below 1 for the other variants. Specifically, we estimate the reproductive number for B.1.1.7 was 1.24 [1.07-1.41] from 01 January until 17 January 2021 and 1.18 [1.06-1.30] from 18 January until 01 March 2021 based on the whole genome sequencing data. From 10 March to 16 March 2021, once B.1.1.7 was dominant, we estimate the reproductive number was 1.14 [1.00-1.26] based on all confirmed cases. For reference, Switzerland applied more non-pharmaceutical interventions to combat SARS-CoV-2 on 18 January 2021 and lifted some measures again on 01 March 2021. CONCLUSION: The observed increase in B.1.1.7 frequency in Switzerland during the study period is as expected based on observations in the UK. In absolute numbers, B.1.1.7 increased exponentially with an estimated doubling time of around 2-3.5 weeks. To monitor the ongoing spread of B.1.1.7, our plots are available online.